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ON  THE  MOTION  OF  NEARLY  SYNCHRONOUS  SATELLITES 

by 

R.  R.  Allan 


SUMMARY 

The  longitudinal  forces  due  to  the  longitude-dependent  part  of  the  Earth's 
gravitational  potential  oan  have  considerable  effects  on  nearly  synchronous 
satellites.  As  a  preliminary  the  special  case  of  the  oiroular  equatorial  orbit 
is  considered.  An  approximate  solution  is  given  for  the  motion  under  the 
combination  of  all  tesseral  harmonics,  and  the  relation  to  the  Jacobi  integral 
is  discussed.  For  general  orbits  the  disturbing  funotion  for  the  combination 
of  all  tesseral  harmonios  is  developed  in  terms  of  the  usual  elliptio  elements 
and  the  resonant  terms  identified  and  isolated.  With  a  suitable  definition  of 
the  mean  longitude,  it  is  shown  that  the  motion  in  lv  -igitude  relative  to  the 
Earth  is  equivalent  to  that  of  a  particle  moving  in  a  one-dimensional  potential, 
provided  only  that  the  eooentrioity  is  small. 
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1  INTRODUCTION 

Although  various  authors  have  disoussed  the  motion  of  nearly  synchronous 
satellites  of  the  earth,  only  Uusen  and  Bailie  and  Morando  have  considered 
other  than  a  oiroular  equatorial  orbit,  the  so-oalled  'geostationary1  satellite. 
The  treatment  given  by  Musen  and  Bailie  is  valid  for  large  eooentrioities  and 
all  inclinations  but  includes  only  the  2  tesseral  harmonio.  Morando,  using 
von  Zeipel's  method  as  modified  by  Hori  to  deal  with  resonanoe,  has  given  the 
stable  positions  and  the  periods  of  libration  for  each  of  the  tesseral  harmonios 
separately  up  to  l  »  4,  m  ■  4.  In  a  previous  paper1-  some  results  were  given  for 
the  oiroular  equatorial  orbit  whioh  hold  for  the  superposition  of  all  tesseral 
harmonios,  and  the  objeot  of  the  present  paper  is  to  extend  these  results  to 
more  general  orbits  and  in  particular  to  nearly  oiroular  but  lnolined  orb  it  s ■ 

2  THE  GRAVITATIONAL  POTENTIAL 

We  will  write  the  gravitational  potential  U(r,0,<p)  due  to  the  earth  at 
distanoe  r  from  the  oentre  of  the  earth,  arid  at  the  point  with  colatitude  6  and 
longitude  9,  in  the  form 

U(r,6,*)  «  UQ(r)  +  U'(r,6,<p)  ,  (l) 

where  UQ  =  p/r  is  the  oentral  potential,  with  p  »  GM  where  M  is  the  mass  of  the 
earth,  and 


U'(r,0,9) 


CO  l 

(w/r)  ^  ^  6)  00. 

vs2  m=1 


.  (2) 


Here  R  is  the  mean  equatorial  radius  of  the  earth,  and  and  9^  are  the 
oonstants  associated  with  the  (6,m)  tesseral  harmonio.  Also  P*(s)  is  the 
associated  Legendre  function  defined  by 


*>■<*)  ■  (1  - ,2>V2  pt(,) 


c"  ll 


9  . 

ds 


(3) 


The  terms  with  m  e  0  would  give  the  sonal  harmonios,  i.e.  the  axially  synuostrio 

part  of  the  field,  and  are  omitted.  Also  the  possible  term  in  (2)  with  l  a  2, 

m  ■  1  must  be  very  small  sinoe  the  axis  of  rotation  of  the  earth  must  very 

5 

nearly  ooinoide  with  the  prinoipal  axis,  and  this  term  is  also  omitted. 
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In  general  the  tea serai  harmonios  lead  only  to  short-period  variations; 
indeed,  values  of  the  oonstants  and  have  been  determined  from  these 
short-period  variations  in  the  orbital  elements  of  olose-earth  satellites.  For 
a  nearly  synohronous  satellite,  however,  the  longitudinal  foroes  due  to  the 
tesseral  harmonios  aot  continuously  in  the  same  sense  and  produoe  long-period 
ohanges  in  the  energy,  the  semi-major  axis,  and  the  mean  motion,  thus  leading 
to  a  libration  in  longitude.  The  main  effeot  cf  the  axially-symmetrio  terms 
whioh  have  been  omitted  from  (l)  and  (2)  (the  only  significant  term  is  the 
seoond  zonal  harmonio  whioh  corresponds  to  the  oblateness  of  the  earth)  is  to 
produoe  a  regression  of  the  orbital  plane  whioh  at  synohronous  height  is  of 
muoh  longer  period  than  the  libration  in  longitude.  The  motion  of  distant 
oiroular  orbits,  inoluding  synohronous  orbits,  under  the  oombinad  effeot  of 
the  earth's  oblateness  and  the  luni-solar  foroes  has  recently  been  dismissed  by 
Allan  and  Cook*’, 

3  CIRCULAR  EQUATORIAL  OBBU 

If  the  orbit  is  oiroular  and  equatorial,  the  satellite  will  appear  to  be 
nearly  stationary  relative  to  the  earth,  apart  from  the  long-period  motion  in 
longitude,  and  this  speoial  oase  is  amenable  to  a  simple  treatment.  We  will 
ignore  the  small  north-south  component  of  foroe  due  to  the  odd  tesseral 
harmonios  (i.e,  those  with  d-m  odd,  whioh  are  antisynmetrio  about  the 
equatorial  plane)  and  assume  that  the  satellite  remains  in  the  equatorial  plane. 
Taking  the  polar  coordinates  of  the  satellite  as  r  and  <p,  and  the  angular 
velooity  of  rotation  of  the  earth  as  nQ,  the  mean  motion  of  the  satellite  is 
n  •  n^  +  f .  Then  the  radial  and  transverse  equations  of  motion  oan  be  written 
as 


9  -  r(no+$)2 


(4) 


£  fr2(n0+  $)!  »  2*(n0+$)+i^  -  F^  ,  (5) 

where  Fr  and  F^  are  the  radial  and  transverse  components  of  the  disturbing 
foroe  derived  from  (2).  We  may  now  expand  about  the  equilibrium  solution  of 
these  equations  for  sero  perturbing  foroe,  namely 


“3 
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2  t 

$  »  0,  p  ■  a0,  n0  ao  *  ^  *  (0 

by  writing  r  ■  aQ  +  p,  and  retaining  only  terms  of  first  order  in  p  and  $. 

This  procedure  leads  to  the  equations 

?  -  3»o  P  “  2a0  no  ♦  "  Pr  *  <7) 

2no  p  +  a0  3  -  •  (8) 

By  eliminating  p  between  (7)  and  (8)  we  oan  derive  the  following  single 
equation  for  9 


[*' 


(It) 


+  n2 
+  n0  <P 


] 


•  3n  P  -  2n 
^09  o 


ft  (V 


dt 


(P_) 

9 


(9) 


Since  d/dt  operating  on  Fy  or  is  equivalent  to  f(d/3r)  +  9(3/09),  the  second 

term  on  the  right  of  (9)  is  of  higher  order  than  the  first  term,  and  may  be 

negleoted.  Similarly  the  third  term  on  the  right  of  (9)  is  smaller  still,  and 

is  also  negleoted.  Also  provided  the  long-period  motion  in  longitude  is  slow 

(iv) 

oompared  with  the  rotation  of  the  earth,  9'  on  the  left  may  also  be  negleoted. 
Finally  this  leaves  the  approximate  equation 

a  B  «  -  3F  .  (10) 

o  <p 

It  is  worth  noting  that  aooording  to  (lO),  the  apparent  longitudinal  accelera¬ 
tion  relative  to  the  earth  in  the  long-period  motion  is  three  times  as  great 
as  the  true  aooeleration  and  in  the  opposite  sense.  Physically  if  the  longitu¬ 
dinal  foroe  is  in  the  direotion  of  motion  of  the  satellite  (F^  >  0),  the  energy 
is  increasing,  so  that  the  semi-major  axis  is  also  increasing  and  the  mean 
motion  is  deoreasing. 

Equation  (lO)  contains  only  the  longitudinal  oomponent  of  the  disturbing 
foroe  in  the  equatorial  plane,  which  is  given  by 


F 


9 


f;  U*(a0,V2,9)  , 


(ID 


where  the  disturbing  potential  is  given  as  a  general  superposition  of  harmonios 
in  (2),  Then  from  (10)  the  equation  of  motion  is 


6 


'*  "  +  3ne  Yj  “  P“(0)  *in  “(<P“  #  (12) 

Since  p“(0)  vanishes  when  (l  -  a)  is  odd,  only  the  even  tesseral  harmonios 
contribute  to  the  longitudinal  disturbing  foroe.  Equation  (12)  nay  be 
integrated  at  onoe  to  give 


1  .2 

2 


y  / 


P?(o)  008  n(<p-  <P,  )  B  constant 
v  cm 


(13) 


Thus  the  long-period  motion  in  longitude  relative  to  the  rotating  earth  is 
equivalent  to  that  of  a  partiole  moving  with  the  velocity  aQ  $  in  the  one¬ 
dimensional  gravitational  potential  -3U'  (ao,x/2,9). 

Deduotlon  from  the  Jaoobi  Integral  -  The  result  (13)  oan  also  be  derived  as  an 
approximate  form  of  the  Jaoobi  integral  of  the  system.  Assuming  the  satellite 
remains  in  the  equatorial  plane  so  that  we  oan  set  6  »  r/2  and  I  ■  0,  the 
Lagranglan  of  the  system  may  be  written  in  terms  of  r  and  9  as 


L  .  £  *2  ♦  j  r2(no  +  $)2  +  \i/r  +  U*(r,V2,«p)  •  00 


Sinoe  the  Lagranglan  is  time-independent,  the  system  possesses  the  Jaoobi 
integral,  H  a  oonstant,  whioh  takes  the  form 


H 


a  “■  f 2  +  ^  r2($2-  n2)  -  v/r  -  U' (r,V2»<p)  *  oonstant 


(15) 


We  will  now  show  that  this  reduoes  to  the  approximate  result  (13)  to  order 
J,  whore  J  stands  for  J^.  First  of  all  from  (13)  $  is  of  order  J^2,  while 
from  (12)  $  is  of  order  J,  and  in  general  the  kth  derivative  of  9  is  of  order 
jV * .  On  the  other  hand  we  oan  write  equation  (k)  in  the  form 


(v*)‘ 


i»  +  0(J) 


.2  .3 


a  + 


0(J) 


(16) 


whioh  leads  to 


r  "  *°  t1  "  3  \  *  61  +  °(j3/2)] 

where  6^  ie  of  order  J«  Consequently 
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(17) 


■3 


7 


P  ■  “  f  a0  */ao1'  *  0( J1/2) }  , 


(18) 


jro  that  p  la  of  order  J^2,  while  p  (whioh  ia  alao  r)  ia  of  order  J  and  p  ia 


of  order  J- 


3/2 


2  2 
In  (15)  £  may  he  negleoted  since  it  ia  of  order  J,  and  on  substituting 

-1 


from  (17)  for  r  and  r  the  remainder  gives 

*  2  no  ao  (1  ’  ’2/no^  f1  "  3  "t  *  9  ^2  +  261  +  °Cj3/2)l 

'.on  J 


2  2 

-  n  a 
0  0 


2 

[1+3  n""51+9  ^2  +  °(J3/,2)J  “  tf’(a0,V2,<p) 


n  constant  +  o(J5/2)  . 

After  multiplying  out  the  first  product  and  adding,  the  terms  involving  6? 
oanoel,  leaving  the  result 

2 

-  n2  a2  ^  ^  +  0(J^2)^  -  U*(ao,V2>«p)  *  oonatant  , 


i«e» 


j  a2  $2  +  3U'(ao»  V2»<p)  *  oonatant  +  0(J^2)  , 


(19) 


whioh  agrees  with  the  previous  result  (13)« 

In  the  light  of  this  consideration  of  orders  of  magnitude,  it  is  olear 
that  (10)  is  derived  using  only  that  part  of  (7)  whioh  is  of  order  J^2,  namely 

-  3n2  p  -  2ao  nQ  $  =  0  +  0(j)  .  (7*) 

In  fact  (7)  is  not  aoourate  to  order  J,  since  among  the  quadratio  terms 
negleoted  the  following  three  are  of  order  Jt 


*o*2  a 


-  2n  p  $ 


+  3n2  p2/afl 


This,  however,  does  not  affeot  the  results  (10),  (12)  and  (13)» 
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4  DEVELOPMENT  OF  THE  DISTURBING  FUNCTION  FOR  GENERAL  ORBITS 

For  a  more  general  orbit,  the  motion  in  longitude  is  still  oontrolled  by 

i. 

the  longitudinal  foroe  experienced  by  the  satellite,  and  it  oan  be  shown  that 
the  motion  must  have  the  same  general  oharaoter  of  an  oscillation  in  longitude 
even  for  orbits  with  considerable  inclination  and  eooentrioity.  However  it  is 
not  easy  to  evaluate  the  mean  longitudinal  foroe  when  the  satellite  is  no 
longer  nearly  stationary  relative  to  the  earth. 

It  seems  simplest  to  return  to  the  disturbing  potential  (2)  and  express 
it  in  terms  of  the  elements  of  the  orbit.  The  shortest,  and  also  the  most 
general,  development  of  the  disturbing  funotion  seems  to  be  that  given  by 
Izsak^.  He  points  out  that  on  rotating  the  coordinate  system  the  dth  degree 
spherioal  harmonio  P°?(ooa  6)  exp(im<p)  must  transform  to  a  linear  combination 

V  f 

of  the  4th  degree  spherioal  harmonios  (oos  6')  exp(im*  9')  where  0'  and  «p* 

are  the  polar  coordinates  in  the  new  system.  If  the  new  ooordinate  system  is 

chosen  with  its  pole  02'  ooinoident  with  the  pole  of  the  satellite  orbit  then 

0'  a  7/2  and  the  assooiated  Legendre  polynomials  P®  (oos  0')  reduce  to  oonstants. 

Moreover  sinoe  P®  (0)  vanishes  when  4  +  m*  is  odd,  the  sum  oontains  only  terms 

with  4  +  m'  even.  The  actual  form  of  the  result  oan  be  derived  from  texts  on 

3  q 

the  applications  of  group  theory  in  quantum  mechanics  (e.g.  Wlgner  or  Rose7), 
end  oan  be  written  as 


4 

p“(oos  0)  exp(im9)  ■  ^  K!^“2^  "PUtC*’  2j)(«  +  f)  ♦  m(0-not))}  ,  (20) 

0 


where 


*4-2  j 
K4m 


,4-n 


LLL&H _  V  2  A  /  2j  \  54-1-2 J 

2*  ji(4-j)t  I.  }  \  k  A«— V 


-2k  _-4+m+2 j+2k 
0 


...  (20 


In  (20),  Q  is  the  right  asoension  of  the  asoending  node,  e»  is  the  argument  of 
perigee,  f  is  the  true  anomaly,  and  nQt  is  the  ourrent  sidereal  time  of 
Greenwioh,  so  that  n-  nQt  is  the  angle  from  the  ourrent  position  of  Greenwich 
to  the  asoending  node  and  <0+  f  is  the  angle  from  the  asoending  node  to  the 
ourrent  position  of  the  satellite  (see  Fig.l).  Also  in  (21),  0  a  oos  1/2  and 
s  a  sin  1/2  where  I  is  the  inolination,  so  that  the  K-  ooeffioienta  depend  only 


3 


9 


on  the  Inclination.  The  form  Kf“2J  quoted  here  la  identloal  to  Isaak* a  k} 

a  Ctt  KM 

exoept  that  the  faotor  ic  haa  been  inoluded  within  the  ooeffioient. 

With  the  aid  of  (20)  the  diaturbing  potential  (2)  oan  be  written  as 


where  the  contribution  of  the  (d,m)th  teaaeral  harmonio  ia 


U 


dm 


(|V*0  Jgffl(Vr)4+1  *  ^ 

...  (23) 


j-0 


Here  #  denotea  the  real  part.  It  has  been  convenient  in  (23)  to  introduoe  the 
quantity  (p*^  given  by 


nt  -  (1 
o 


M 


whioh  is  the  ourrent  longitude,  relative  to  the  ascending  node,  of  the  meridian 
assooiated  with  the  (d,m)th  tesseral  harmonio. 

The  expression  (23)  still  oontains  the  true  anomaly  and  must  be  developed 
in  terms  of  the  mean  anomaly  M  ao cording  to 

00 

(r/a)""6"1  exp(ipf)  *  x-d-1,p  exp(iqM)  ,  (25) 

qc-00 

which  introduoes  the  Hansen  ooeffioients  These  ooeffioients  depend  only 

on  the  eooentrioity  and  oan  be  expanded  as  power  series  in  e,  and  we  will  re o all 
the  following  properties: 

(i)  If  p  ■  q,  o  1  +  0(e2), 

(ii)  If  p  *  qf  the  leading  term  in  X®*^  is  of  order  |p-  q|  in  e. 

Then  from  (23)  and  (25) 

®{.  •  <«*>  a  Y,  Z  rU>  x;*-’'**21  «p[i[U-2j)  .♦  , 

i  q 


•  •• 


(26) 
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whioh  oompletes  the  development  of  the  disturbing  funotion  in  terms  of  the 
elements  a,e,I,n,w  and  x'*  where  x*  i*  the  modified  mean  anomaly  at  the  epoch, 
defined  by 


dt  +  x' 


To  piok  out  the  terms  whioh  will  be  important  for  synchronous  orbits,  we 
note  that  the  rate  of  ohange  of  the  argument  in  (26)  is 

(l-  2j)  &  +  q(n+  x')  -  m(nQ-  ft)  «  qn  -  m  nQ  .  (28) 

Sinoe  n  “  nQ,  the  resonant  terms  are  those  with  q  ■  m,  and  the  resonant  part  of 


the  potential  oan  be  written  as 


U*B  =  (hAO  JgB(Va)4+1  ft  Yt  4;2d  e*p[i[(*-2j)  «  +  • 

i 

...  (29) 

Sinoe  the  Hansen  ooeffioients  in  (29)  are  of  order  j-6-  m-  2j|  in  e,  the  lowest 
order  of  term  whioh  aotually  ooours  depends  on  whether  (d-m)  is  even  or  odd: 

(i)  If  (l-m)  is  even,  the  lowest  order  term  is  given  by  2J  a  g-  m  and 
is  independent  of  e. 

(ii)  If  (l-m)  is  odd,  the  lowest  order  terms  are  those  of  order  e 
given  by  2 j  a  l  -  m  ±1 . 


I  (»)  ^  gFq  Z*L>)  j  *3  lOg  J  go  l) 


It  seems  desirable  to  oonsider  in  more  detail  the  term  oontalning  Jg  g 
whioh  oorresponds  to  the  elliptioity  of  the  earth's  equator.  Although  the 
ooeffioients  are  not  very  well  determined  as  yet  it  is  probable  that  2 
is  the  largest  coefficient,  and  fairly  oertain  that  it  will  give  the  dominant 
term  at  synchronous  height  sinoe  U^B  decreases  with  distanoe  as  r"^+1\  Thu: 
from  (29),  evaluating  the  K-  ooeffioients  aooording  to  (21 )  and  expanding  the 
Hansen  ooeffioients  up  to  order  e\  the  resonant  part  of  the  disturbing  funo¬ 
tion  for  this  harmonio  takes  the  form 
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u$#2  "  3(wA0  J2|2(V*)3 

x[(l  -  |  e2  +  o*  ooa  2(ll- ^>2  +  «) 

+  (f  ®2  +  2  +  *•*)  °2  ®2  008  2(M“  **,2^ 

+  (§£  +  ...^  a^  008  2(M-<p*^-0)Jj  .  (30) 

This  result  oould  also  be  written  down  from  the  work  of  Cook1®  who  has  given  an 
explioit  development  of  the  disturbing  function  up  to  l  =  m  =  4.  He  was*  how¬ 
ever*  only  concerned  with  the  'long-period'  terms  for  non-resonant  conditions 
and  eventually  discarded  the  terms  whioh  are  important  here. 

The  rates  of  ohange  of  all  the  elements*a*e*I*fl,a>  and  x'  can  now  be  derived 
from  (30)  by  Lagrange's  planetary  equations11.  As  will  be  shown*  the  motion  of 
synchronous  satellites  is  determined  almost  entirely  by  the  variation  in  the 
semi-major  axis,  and  the  variations  of  all  the  other  elements  are  unimportant. 
Consequently  we  give  to  order  e^  only  the  equation  for  d«/dt.  The  Lagrange 
equation  for  a*  namely 

d*/dt  -  (2/na)(dU'/dx')  *  (31) 

leads  to  the  result 


Clearly  the  most  useful  variable  to  desoribe  the  notion  is  simply  the 
oonbination 


M  -  fg  2  +  o>  »  (M+co+n)  -  nQ  t  -  2 
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whioh  ooours  in  the  argument  of  the  resonant  term  of  sero  order  in  e  in  (30). 

The  slowly-varying  quantity  (u+  o>+n)  -  nQt  defines  the  mean  longitude  of  the 
satellite  relative  to  the  earth,  and  we  will  write 

n  dt  -  nQ  t  +  s*  .  (33) 

where  e'  ■  x'  +  »  +  0  is  the  modified  mean  longitude  at  the  epooh. 

It  is  fairly  easy  to  interpret  <p  in  terms  of  the  positions  of  the  nodes.  For 
an  exaotly  synchronous  satellite  (n  =  nQ)  the  ground-traok  on  the  earth  is  a 
fixed  ourve  whioh  degenerates  to  a  symmetrical  figure-of-eight  centred  on  the 
equator  if  the  orbit  is  also  oiroular  (see  Fig.2).  For  an  eooentrio  orbit  the 
true  anomaly  is  -u  at  the  asoending  node  (t  a  tN),  and  (x-o>)  at  the  descend¬ 
ing  node  (t  a  tN, ),  so  that  from  Kepler's  equation  the  mean  anomaly  takes  the 
values  (—  <0-4-0  sin  o>)  and  (x-  w-  e  sin  00)  at  the  asoending  and  descending  nodes 
respectively.  Using  these  values  in  (33)  the  longitudes  of  the  asoending  and 
descending  nodes  relative  to  the  earth  are  respectively 

D  -  nQ  tN  »  9  -  e  sin  w  , 

(fi+x)  -  nQ  tN,  a  fie  sin  co  , 

i.e.  5  is  given  by  the  mean  of  the  asoending  and  descending  nodes  on  the  earth. 
In  praotioe  the  satellite  is  not  exaotly  synchronous  so  that  9  is  slowly- 
varying  and  the  ground-traok  oh&nges,  but  the  value  of  9  at  any  instant  oan 
still  be  determined  by  Interpolation 

Using  the  mean  longitude  9  the  rates  of  ohange  of  all  the  elements  for  a 
nearly  oiroular  orbit  are  as  follows: 

dv'dt  a  -  12  na  J2^2(R/*)2  o*  sin  2(9-  9^)  +  0(e2) 

de/dt  a  3ne  J^2(B/*)2  sin  2(9-  92^2)  -  9  o2  s2  sin  2(9-  92^-«#)|  +  0(e5) 

a  6n  J2^2(R/*)2  0^  s  sin  2(9-  9^)  +  0(e2) 

diVdt  a  -  3n  J2^2(R/»)2  o2  oos  2(9-  92^2)  +  0(e2) 


(ll+to+n)  -  n 
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d«/dt  ■  3n  J2>2(V*)2  f"  o2(1  +  3  o2)  oos  2 ($-l>2>2) 

+  9  o2  s2  oos  2(5- <p2^2-»)}  +  0(o2) 

dx'/dt  a  3n  J2j2(V*)2  [11  00*  2(5-  <f2^2) 

-  9  o2  s2  oos  2(5-  *2,2“  +  °(®2)  *  (34) 

There  are  terns  in  the  last  two  equations  of  (34)  whioh  are  of  aero  order  in  e 
yet  still  depend  on  o>,  but  these  terms  oanoel  on  forming  the  equation  for  5* 

Thus  from  (33) 

d5/dt  a  (n-nQ)  +  de'/dt  , 

and  differentiating  onoe  mare, 

d25/dt2  a  di/dt  +  d2e'/dt2  .  (35) 

Then  troa  the  first  of  (34), 

dVdt  »  +  I8n2  J2^2(B/»)2  oV  sin  2(5-  ,^2)  ♦  0(e2)  .  (36) 

Likewise  from  the  last  three  of  (3(4), 

de'/dt  ■  6n  J^2(b/*)2  o2(l>  o2-l)  oos  2(5-  +  °(®2)  »  (37) 

and  differentiating  onoe  more 

d2s*/dt2  •  -  12n  d^dt  J2^2(V«)2  o2(4  o2-l)  x  sin  2(*-f2>2)  *  °(®2) 

...  (38) 

Sinoe  d$/dt  «n,  the  second  term  on  the  right  of  (35)  is  negligible  oompared 

with  the  first  term.  In  faot  (36)  is  of  order  J,  while  (38)  is  of  order  J^2 

sinoe  dq/dt  is  of  order  J1/2.  To  lowest  order  the  notion  in  mean  longitude  is 

controlled  only  by  the  variation  in  the  semi-major  axis,  and  if  the  seoood  term 

2 

in  (35)  is  negleoted  the  equation  of  motion  oan  be  written,  to  within  order  a  , 


as 
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d2^at2  ■  2  ^ 

where 

k2  -  36n2  J2<2(R/a)2  o4  .  (40) 

There  are  olearly  four  positiona  of  equilibrium;  the  two  positions 
<p  *  <P2  2  ±  -r/2  (on  the  minor  axis  of  the  earth's  equatorial  section)  are 
stable,  while  those  at  9  a  <P2  2  4114  i  a  <P2  2  +  *  (on  the  maJor  axi®)  are 
unstable*  To  integrate  (39)  it  is  more  convenient  to  ohange  the  origin  to 
one  or  other  of  the  two  stable  positions  by  writing 

♦  ■  5  -  *>2,2  *  ^2  9 


so  that  (39)  becomes 

’♦  *  -  ^  k2  sin  2  ijr  •  (41 ) 

The  first  integral  of  (41 )  oan  be  written  as 

$r2  -  k2  oos2  ijr  *  oonstant  •  $2  -  k2  ooa2  ,  (42) 

where  \jrQ  and  are  the  initial  values*  The  satellite  will  be  captured  in  an 
effeotive  potential  well  and  will  oscillate  about  one  or  other  of  the  two 
points  of  stable  equilibrium  provided  that 

♦0  <  00*2  ^0  *  (43) 

If  this  oondition  is  satisfied,  ijr  will  osoillate  between  the  limits  ±  ^  say, 
where 


k2  sin2  t a  ■  ♦  42  sin2  ,  (l»4) 

and  (42)  oan  be  written  as 

<r2  a  k2(sin2  -  sin2  $)  •  (45) 


From  (45)  the  period  of  a  complete  oeoillation  is 


/3 
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*n  V2 

T  ■  | *  I  (sin2  $B-sin2  m  k  J  (1  -  Bin2  f/m  sin2  u)"1^2  4u  * 

o‘  oJ 

•••  (46) 

The  result  can  be  written  in  the  alternative  forms 


T  a  £K(sin*B)  =  2f  ^(^jljsin2^) 


(47) 


where  K(sin  \|/m)  la  a  oomplete  elliptic  integral  of  the  first  kind,  and  ^  is 

the  ordinary  hypergeometrio  funotion.  If  the  amplitude  of  the  oscillation  is 

small,  the  period  is  approximately  2 x/k,  which  is  about  780  days  for  an  orbit 

of  low  inclination  if  we  take  J,  ,  «  2  x  10*  .  The  period  increases  with  the 

amplitude  of  the  osoillation  and  approaches  infinity  as  *  approaches  r/2, 

m 

To  return  to  the  variation  of  the  remaining  elements  in  (34),  it  is  dear 

that  de/dt  is  completely  negligible.  Both  dl/dt  and  d^/dt  must  certainly  be  less 
2 

than  3n  J„  9(u/a)  whioh  is  approximately  5  x  lO-3  degrees/day  taking 
J2  2  “  2  x  10  •  The  inclination  oscillates  as  the  satellite  performs  a  libra- 

tion  around  a  position  of  stable  equilibrium.  The  remaining  quantities  do/dt 
and  dx'/dt  are  rather  larger,  and  also  more  complicated;  in  particular  there  are 
parts  of  sero  order  in  e  whioh  depend  on  the  argument  of  perigee.  However  for 
the  more  meaningful  combination  e 1  given  in  (37), 


de'/dt  <  I8n  J2>2(h/»)2  , 

whioh  is  approximately  3  *  10*^  degree/day,  again  taking  J,  ,  *  2  x  10**^.  Far 
the  first  tern  on  the  right  of  (35)#  we  have  from  (36)  that 

dn/dt  <  18n2  J-  „(B/a)2  , 

whioh  is  approximately  2  x  10  •*  degree/day  .  From  (36)  and  (38)  the  ratio 
between  the  two  terms  on  the  right  of  (35)»  taking  I  a  0  for  simplicity,  is 
(2/n)  dq/dt  +  0(e2).  Now  from  (45)  d<j/dt  a  k,  where  k  «  10**^  seo"*^  for 
J2,2  “  2  x  10”  »  80  th*t  the  ratio  is  at  most  about  2.5  x  10”^,  and  the  seoond 
term  on  the  right  of  (35)  i a  only  a  small  oorreotion. 
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6  MOTION  UNDER.  THE  COMBINATION  0?  ALL  TESSERA!  HARMONICS.  FOR  e  SHALL 

Prom  the  development  of  Seotion  4  end  in  particular  from  (29)  at  seq,  the 
resonant  part  of  the  oomplete  disturbing  funotion  to  lowest  (sero)  order  in  e 
may  be  written  as 

Ei 

J4m(R/a)4+1  kJb  oos  m(M+u>-  <p$ffl)  +  0(e)  .  (48) 

Here  the  dash  denotes  that  the  summation  is  only  over  the  even  tesseral  harmonios 
the  part  of  order  e  arises  from  the  odd  tesseral  harmonios.  In  writing  (48),  we 
have  used  the  properties  of  the  Kansan  ooeffioients  and  also  the  faot  that, 
aooording  to  (21 ),  the  K-  ooeffioients  are  real  for  4-  m  even.  This  disturbing 
funotion  is  then  adequate  to  describe  the  motion  for  orbits  at  any  inclination 
provided  the  eooentrioity  is  small. 

From  (20)  it  is  olear  that  when  I  *  0,  0  *  r/2 ,  and  consequently 

m  for  I  ■  0.  Thus  it  is  possible  to  write,  when  4-  m  is 

even, 

-  Dt.(I)  p!(0)  •  <w> 

where  each  ■  1  when  1*0.  Also  from  (24)  and  (33) 

M  +  ®  t  (50) 

so  that  the  resonant  part  of  the  die  -urbing  funotion  takes  the  form 

N*  -  (mA)  Y_  J4m(R/a>4+1  P4(0)  D4m(l)  008  +  »  (51) 

where  the  sum  is  over  the  even  tesseral  harmonios.  The  ooeffioients  D^m  for 
the  five  even  tesseral  harmonios  up  to  4  *  m  *  4  are  easily  found  from  (21)  or 
from  the  explicit  development  of  Cook10  and  are  given  in  Table  1  together  with 
the  values  of  P*(0). 
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TABLE  1 

The  ooeffioients  D^a  up  to  t  ■  m  ■  4 


Harmonios 

*J<°) 

(2,2) 

3 

V 

-  0^ 

(3,0 

-3/2 

D3.1 

>1  o2  (1  -  10  s2+  15 

(3,3) 

15 

D3.3 

6 

3  0 

(4,2) 

-15/2 

\,2 

-  o^  (1  -  14  s2  +  28 

(4,4) 

105 

-  o8 

As  before  the  equation  of  motion  for  $  is  given  by  (35)  where  the  first 
and  seoond  terms  on  the  right  are  of  order  J  and  J^2  respectively.  Negleoting 
the  higher  order  term  and  oombining  (35)  and  (31 ), 

d2f/dt2  «  dn/dt  »  -  (3/a2)  dU*/d9  (52) 

where  3/3 x'  in  (31 )  has  been  replaced  by  3/dlp  on  account  of  (50).  Then  from 
(51)  and  (52),  * 

_  t 

d25/dt2  »  3n2  )  »  ^(H/a)4  Pj(0)  Dto(l)  sin  »(9-^#)  +  0(e) 

...  (53) 

In  finding  the  equation  of  motion  for  the  mean  longitude,  we  have  taoitly 
assumed  that  the  odd  tesseral  harnonios,  whioh  are  lumped  together  as  the  terms 
of  order  e  in  the  potential,  only  contribute  terms  of  order  e  to  (53).  This  is 
fairly  clear  physioally  since  the  odd  harmonios  are  antisymmetric  about  the 
equator*  However  we  ought  to  oonsider  the  effoot  of  the  odd  harmonios  on  all 
orbital  elsments  as  we  have  done  in  (34)  for  the  Jg  g  term.  Ooing  back  to  (29) 
and  using  Lagrange's  planetary  equations,  it  is  easily  shown  that  the  odd 
harmonios  contribute  terms  of  order  e  to  de/dt,  dl/dt  and  diVdt*  As  is  to  be 
expeoted  the  odd  harmonios  are  acre  effective  in  o hanging  the  position  of 
perigee  and  give  terms  of  sero  order  in  e  in  de/dt  and  of  order  e  ^  in  da/dt  and 
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TABLES  1 

The  eoefflolenta  up  to  l  »  a  ■  4 


Harmonios 

*JCo) 

(2,2) 

3 

D2,2  * 

o* 

(3,0 

-3/2 

“3,1  * 

o2  (1-10  s2+  15  S4) 

(3,3) 

15 

”3,3  ’ 

6 

0 

(4,2) 

-15/2 

v  • 

oK  (1  -  14  s2  +  28  Bk) 

(4,4) 

105 

V  * 

o8 

As  before  the  equation  of  motion  for  <p  is  given  by  (35)  where  the  first 
and  seoond  terms  on  the  right  are  of  order  J  and  J^2  respectively.  Negleoting 
the  higher  order  term  and  combining  (35)  end  (31), 

d2y/dt2  «  dn/dt  »  -  (3/a2)  dU*/*¥  (52) 

where  d/dx'  in  (31 )  has  been  replaced  by  d/d$  on  account  of  (50).  Then  from 
(51)  and  (52),  • 

d2<j/dt2  »  3n2  y  a  0^(100*  pJ(O)  0^(1)  sin  ■(»-  ^#)  +  0(e) 

...  (53) 

In  finding  the  equation  of  motion  for  the  mean  longitude,  we  have  taoitly 
assumed  that  the  odd  tesseral  harmonios,  whioh  are  lumped  together  as  the  terms 
of  order  e  in  the  potential,  only  contribute  terms  of  order  e  to  (53).  This  la 
fairly  dear  physioally  since  the  odd  harmonios  are  antisymmetric  about  the 
equator.  However  we  ought  to  oonsider  the  effeot  of  the  odd  harmonios  on  all 
orbital  elsments  as  we  have  done  in  (3 4)  for  the  Jg  ^  i®rB»  Going  baok  to  (29) 
and  using  Lagrange's  planetary  equations,  it  is  easily  shown  that  the  odd 
harmonios  contribute  terms  of  order  e  to  d</dt,  dl/dt  and  d£}/dt.  As  is  to  be 
expeoted  the  odd  harmonios  are  more  effeotive  in  0 hanging  the  position  of 
perigee  and  give  terms  of  aero  order  in  e  in  de/dt  and  of  order  e”^  in  dtq/dt  and 
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dx'/dt  separately.  However  the  contribution  to  the  nore  meaningful  combination 
de'/dt  is  easily  shown  to  be  of  order  e.  Thus  the  odd  teaseral  harmonios 
contribute  only  terms  of  order  e  to  (53).  On  the  other  hand  all  the  even 
tesseral  harmonios  give  terms  in  da/dt,  etc  of  the  same  form  as  in  (3k)  • 

The  result  given  in  (53)  is  similar  in  form  to  (12)  exoept  that  the 

longitude  <p  in  (l 2 )  has  been  replaoed  by  the  more  general  quantity  <p»  and  (53) 

also  oontains  the  faotors  D,  (l).  As  for  the  geostationary  oase  the  motion  in 

cm 

9  is  equivalent  to  that  of  a  particle  moving  in  a  one-dimensional  potential  and 
in  prinoiple  the  satellite  could  be  oaptured  in  any  potential  minimum.  Of 
oourse  from  the  estimates  which  have  so  far  been  obtained  for  the  ooeffioients 
it  is  very  probable  that  the  potential  is  dominated  by  the  (2f2)  harmonio 
whioh  gives  just  two  minima  in  the  potential.  However  it  is  just  possible  that 
the  superposition  of  the  other  tesseral  harmonios  oould  give  rise  to  further 
minima,  particularly  by  splitting  the  prinoipal  minima  of  the  (2,2)  harmonio 
into  a  number  of  separate  minima. 

If  the  ooeffioients  DgB(l)  mre  assumed  to  be  sufficiently  constant  over 
the  period  in  question,  (53)  could  be  integrated  to  give 

Et 

J4m(V»0)4  p*(°)  Dga(l)  008  ■(?'  Tto)  *  constant  +  0(e) 

...  (5k) 

However  this  integral  is  more  suspect  and  less  useful  than  the  corresponding 
result  (13)  for  a  nearly  equatorial  satellite,  sinoe  the  inclination  of  the 
orbit  is  changing  with  time  •  It  is  true  that  the  inclination  of  an  initially 
equatorial  orbit  also  ohanges  with  time,  but  dD^y'dl  vanishes  for  an  equatorial 
orbit  although  not  for  an  inolined  orbit.  In  faot  two  of  the  ooeffioients  in 

Table  1  change  considerably  with  inclination;  D_  .  becomes  negative  when  the 

0  'f1  0 

inclination  exoeeds  1*0  and  beoomes  positive  again  above  95  .  Likewise  „ 

has  zeroes  at  3k°  and  80°  inclination.  The  remaining  ooeffioients  in  Table  1 

are  positive  for  all  inclinations. 

Although  the  potential  is  very  probably  dominated  by  the  (2,2)  harmonio, 
it  is  still  interesting  to  consider  the  motion  under  the  general  (d,m)  tesseral 
harmonio  alone.  From  (53),  if  only  the  (d,m)  harmonio  is  present,  there  will 
be  2m  positions  of  equilibrium  divided  into  two  sets  each  of  m  as  follows; 

(i)  f  -  9lm  +  2p 

(ii)  ♦  ■  fga  ♦  (2P+ 0 
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where  p  ■  0,1 ,2,...(m- 1 ).  If  the  ooeffioient  P®(0)  D^o(l)  is  positive, 
the  former  are  unstable  and  the  latter  are  stable,  and  vioe-versa  if  the 
ooeffioient  is  negative.  In  any  oaae  the  ooeffioient  oan  be  made  positive  by 
suitably  redefining  <pga  so  as  to  ohange  the  sign  of  Assuming  the 

ooeffioient  is  positive,  the  substitution 

5  -  <p£m  »  2t/m  +  (2p+l)  »  (55) 

ohanges  the  origin  to  a  stable  point,  and  reduoes  the  equation  of  motion  to 

5  =  -  £  k*#  sin  2*  +  0(e)  ,  (#) 

where 

4,  -  pPo)  •  (57) 

Since  (56)  is  equivalent  to  (4l),  the  equations  (A2)  to  (V7)  hold  with  k 
replaoed  by  k^ffl. 

7  CONaUSIONS 

As  is  now  well-known  the  very  small  longitudinal  forces  due  to  the  slight 
asymmetry  of  the  Earth  about  its  axis  oan  have  considerable  effeots  on  nearly 
synchronous  satellites.  Beoause  these  foroes  oan  aot  continuously  in  one  sense, 
the  energy  of  the  satellite  in  its  motion  around  the  Earth  slowly  ohanges.  As 
a  preliminary  some  further  consideration  has  been  given  in  the  first  part  of 
this  Report  to  the  speoial  oe3e  of  a  satellite  in  a  oircular  equatorial  orbit. 
Here,  as  has  been  shown  before,  the  motion  of  the  satellite  relative  to  the 
Earth  is  simply  that  of  a  particle  moving  in  a  one-dimensional  potential  which, 
apart  from  a  constant  faotor,  is  the  gravitational  potential  at  synchronous 
height  on  the  equatorial  plane.  This  involves  only  the  even  tesseral  harmonics 
(i.e.  those  with  l-  a  even)  since  the  odd  harmonics  vanish  on  the  equatorial 
plane.  Almost  oertainly  the  principal  contribution  is  from  the  term  involving 
Jg  g  whioh  corresponds  to  an  elliptioity  of  the  Earth's  equator.  If  the 
satellite  is  sufficiently  near  to  synchronism  it  oan  be  oaptured  and  osoillate 
within  a  trough  of  this  potential. 

In  the  main  part  of  this  Report  these  results  are  extended  to  orbits  at 
general  inclination  and  which  are  not  neoessarily  oircular.  In  its  motion 
relative  to  the  Earth  such  a  satellite  performs  a  figure-of-eight,  whioh  is 
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symmetric  about  the  equator  if  the  orbit  is  oiroular  but  is  otherwise  a 
distorted  three-dimensional  curve.  Tne  mathematical  method  adopted  is  to 
expand  the  disturbing  funotion  in  terms  of  the  usual  elliptio  elements  and  to 
retain  only  the  'resonant'  or  long-period  terms.  Special  attention  has  been 
given  to  motion  under  the  Jg  2  term  alone*  Finally  the  analysis  is  applied  to 
motion  under  the  combination  of  all  tesseral  harmonios  provided  only  that  the 
eooentrioity  is  small.  It  is  found  that  the  variation  of  the  mean  longitude 
relative  to  the  Earth  is  a  simple  extension  of  the  corresponding  result  for 
the  oiroular  equatorial  orbit,  and  involves  only  the  introduction  of  an  extra 
faotor  D^ffl(l),  depending  only  on  the  inclination  I,  for  each  tesseral  harmonio. 
As  before,  the  motion  in  mean  longitude  is  equivalent  to  that  of  a  partiole 
moving  in  the  appropriate  one-dimensional  potential,  and  the  coefficients 
D^l)  express  how  the  satellite  in  its  figure-of-eight  motion  samples  the 
longitudinal  force  due  to  eaoh  of  the  tesseral  harmonios.  Only  the  even 
tesseral  harmonios  are  Involved  since  the  odd  harmonios  are  anti symmetric  about 
the  equator  and  their  effeots  oanoel  out. 

The  analysis  presented  here  should  permit  the  determination  of  the 
ooeffioients  of  the  even  tesseral  harmonios  from  observations  on  synchronous 
satellites.  For  at  least  some  of  these  ooeffioients,  the  aoouraoy  should  be 
higher  than  oan  ourrently  be  obtained  from  olose-orbit  satellites.  Apart  from 
the  oonverse  application  in  predioting  or  controlling  the  motion  of  future 
synchronous  satellites,  such  results,  possibly  taken  in  conjunction  with  olose- 
earth  results,  should  improve  our  knowledge  of  the  shape  of  the  geold  and  the 
Earth's  external  gravitational  field. 


The  author  is  most  grateful  to  D.  G.  King-Hele  and  G.  E.  Cook,  who  read 
this  paper  in  manuscript  form,  for  their  helpful  suggestions. 
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Fig.l 


FIG.  I  PROJECTION  OF  THE  SATELLITE  ORBIT 

ON  THE  UNIT  SPHERE 
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SPA/P/1335. 


FIG  2  GROUND  TRACK  OF  SYNCHRONOUS  SATELLITE  ON 

THE  ROTATING  EARTH. 
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